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THERMODYNAMIC  EXTRAPOLATION  OF  ROCKET  ENGINE  PERFORMANCE 
PARAMETERS 


Ge  Minglong 


Abstract 


In  this  paper,  extrapolation  formulas  of  rocket 
engine  performance  parameters  are  established  utiliz¬ 
ing  one  isentropic  reference  line  and  two  partial 
derivatives  for  fixed  enthalpy.  These  formulas  can 
be  used  for  the  extrapolation  calculations  of  some 
parameters  such  as  specific  thrust,  characteristic 
velicity,  nozzle  area  ratio,  and  the  thermodynamic 
parameters  at  the  combustion  chamber  and  the  nozzle 
exit  when  changes  occur  in  the  initial  enthalpy  of 
the  propellant,  combustion  chamber  pressure  and  nozzle 
exit  pressure.  The  formulas  in  this  paper  are  simple, 
convenient  to  use  and  have  relatively  high  accuracy. 


I.  INTRODUCTION 


Methods  of  thermodynamic  calculation  of  the  theoretical 
performance  of  rocket  engines  are  classic.  These  are  introduced 
in  reference  1  and  other  literature  dealing  with  rocket  engines. 
These  diverse  calculations  are  very  useful  for  rocket  engine 
design,  but  consume  quite  a  lot  of  machine  time.  It  is  also 
inconvenient  to  put  a  large  quantity  of  data  into  book  form. 
Therefore,  the  problem  of  a  calculation  method  for  thermodynamic 
extrapolation  of  theoretical  performance  parameters  of  rocket 
engines  is  one  that  has  been  raised  for  a  long  time.  Formulas 
for  the  calculation  of  thermodynamic  extrapolation  have  been 
presented  [2-5] . 


In  recent  years,  with  the  popularity  of  pocket  calculat¬ 
ors,  it  has  become  more  convenient  to  use  extrapolation  to 
determine  the  theoretical  performance  of  rocket  engines.  This 
paper,  therefore,  proposes  a  more  accurate  extrapolation  formula 
based  on  existing  extrapolation  formulas. 


II.  ISENTROPIC  REFERENCE  LINES 


For  the  mixing  ration  0/F  is  any  propellant  composition 

with  a  definite  value.  Because  its  element  composition  is 

fixed,  no  matter  what  the  numerical  values  (H  or  H  )  of  the 

coo 

initial  enthalpy  of  the  propellant  affected  by  such  factors  as 

ambient  temperature  and  regardless  of  adopting  a  combustion 

chamber  pressure  (p  ,p  ,  or  P  )  which  is  too  high, and  a  nozzle 

exit  pressure  (p  ,p  . ,p  ,,  or  p)  which  is  too  high,  the  ideal 

wO  c  c  c 

thermodynamic  process  of  a  rocket  engine  -  the  isentropic 

limiting  equilibrium  expansion  process,  can  both  be  shown  as 

different  isentropic  lines  c  t  e  e“,cte,...  etc.  on  the  same 

o  o  o 

enthalpy-entropy  diagram,  as  shown  in  fig.  1. 

For  the  initial  enthalpy  value  of  given  propellants  (H  ) 

and  combustion  chamber  pressure  (p  ) ,  the  thermodynamic  para- 

CO 

meters  of  the  combustion  chamber  state  point  c  and  the  throat 

o 

state  point  tQ  can  be  determined  through  accurate  thermodynamic 
calculation.  Several  nozzle  exit  pressures  are  then  selected 
and  the  thermodynamic  parameters  of  these  nozzle  exit  state 
points  (co,c",...)  are  calculated  on  the  basis  of  isentropic 
conditions. 


On  isentropic  line  c  tec",  there  are  several  points  with 

O  O  O 


known  values  for  the  parameters  of  pressure  ratio  p  /p,  pressure 

c 


p,  temperature  T,  enthalpy  H,  molecular  weight  M,  constant  pres¬ 


sure  specific  heat  c^  and  the  values  for  the  partial  derivatives 
'  (see  references  1  and  4  for  the  computation 
formulas  and  "method  of  these  two  partial  derivatives) .  This 


isentropic  line  can  therefore  become  the  reference  line  for 
extrapolation  calculations.  The  parameters  of  a  typical  isen¬ 
tropic  reference  line  are  given  in  Table  1.  The  parameters 
given  in  the  top  half  of  this  table  are  taken  from  reference  5. 


On  an  isentropic  reference  line,  apart  from  the  known  para¬ 
metric  points,  the  thermodynamic  parameters  of  each  point  can 


approximately  be  calculated  by  using  the  following  interpola¬ 
tion  formulas: 

F  -  Qlap  1 

(AF).  -  ?(%),  J  _  U) 

In  the  formulas,  F  expresses  the  parameters  T,H,M  etc.,  Q  is  the 
interpolation  coefficient  determined  on  the  basis  of  the  para¬ 
meters  of  the  two  adjacent  known  points. 

For  the  isentropic  reference  line  in  Table  1,  the  results  of 

the  calculation  of  the  interpolation  coefficients  (^,Q  and  Q 

a  M  H 

are  given  in  the  lower  half  of  the  table. 

III.  ISENTHALPIC  PARTIAL  DERIVATIVES  AND  ISENTHALPIC  RELATIONAL 
EXPRESSIONS 

By  referring  to  related  formulas  given  in  references  4  and 
5,  the  following  formulas  are  used  to  express  the  two  isenthal- 
pic  partial  derivatives  in  this  paper: 


9ln?/s  c,M\d\nT/r 

D  -  _  -RfdlnM  V  (dlnM  \ 

_  \  dinp  / h  cpM\  91 aT  >*  \  9ln p  /  r 


In  the  formulas,  R  is  a  general  gas  constant  (1.98726  Kcal/Kg 
mole  K) . 

For  the  isentropic  reference  line  in  Table  1,  the  values  of 
D7'DM'^DT  and  ®DM  which  were  calculated  by  using  the  above  ex¬ 
pressions  are  listed  in  the  lower  half  of  the  table. 


The  relational  expressions  of  temperature  and  molecular 
weight  on  the  isenthalpic  line  are 

-  Dr( Atop)*  (♦) 

,  (AlnM),  -  Dm(*l ap)B  (5) 


By  substituting  the  expression  of  the  first  law  of  thermo¬ 
dynamics  and  the  state  equation  in  the  differential  expression 
of  entropy,  we  obtain 

-v  -  r  at 

For  the  isenthalpic  line  (d!!=0)  ,  from  this  expression  we  obtain 


IV.  THE  ESTABLISHMENT  OF  AN  EXTRAPOLATION  FORMULA 


1.  Combustion  Chamber  Thermodynamic  Parameters  and  Character¬ 
istic  Velocity  Formula 

When  the  propellant  initial  enthalpy  H^fcH  and  combustion 

chamber  pressure  pc*pco,  the  parameters  for  c'  in  fig.  1  can  be 

calculated  from  formula  (1)  based  on  known  data  from  ^H  =H  -H 

c  c  co 

and  point  cq.  The  parameters  of  point  c  can  then  be  calculated 
from  formulas  (4)  and  (5)  based  on  the  parameters  of  point  c' 
and  p=pc  . 

As  a  result,  the  formulas  for  combustion  chamber  tempera¬ 
ture  and  molecular  weight  are  obtained 


*  ;.S7!  .>  >nM9i  »**  •  r  -  ,vQ9tef\w:i  fed 

i  -Si 


*  <1.  r  ■ 


Because  the  characteristic  velocity  c*  forms  a  direct  ratio 
with  <^TC/M~,  the  extrapolation  formula  of  the  characteristic 
velocity  can  be  written  as 


*9  Was 


2.  Formula  of  Nozzle  Exit  (p=p  )  Thermodynamic  Parameters 

c 


Formula  (6)  is  applied  between  c  and  c'  and  between  e  and 


e 


t  . 


lafe—Muke* 

ft  Mt>  pt 


By  substituting  formula  (1)  into  this  expression,  we  obtain 

Qh  o..  / 


In  ..  In  _£*.  + 


(10) 


Qh 


In  the  formula,  M  is  the  approximate  value  of  M  ,.  It  is  the 

CO  c 


known  molecular  weight  of  point  PssPco  near  pc ,  on  the  isentropic 

reference  line.  The  initial  value  of  p  ,  can  be  Calculated 

c 


using  the  following  expression  in  order  to  select  pcQ  so  as  to 


determine  M 


co 


OO 


Q£  -  Pt,  ■ 

i  _  .  • 

Key:  (1)  Initial. 

This  expression  is  approximately  obtained  from  formula  (10) , 


where  Mc„  is  the  known  molecular  weight  of  the  lowest  point  of 


the  isentropic  reference  line. 


The  interpolation  coefficients  QT  ,QdT,Qm,Qdm  and 

formulas  (7) ,  (8) ,  (10)  and  (11)  use  values  in  the  vicinity  of 

p  .  When  /lH  ^  0 ,  the  mean  value  of  c  -t  is  used  and  when 
*co  c  '  o  o 

A  H  >0  ,  the  value  of  p  >p  is  used, 
c  co 

On  the  basis  of  the  known  data  of  e  and  the  value  of 

o 

lnpc 1  calculated  from  for. aula  (10) ,  the  parameters  of  point  e* 

Pco 

can  Be  calculated  from  formula  (1) .  On  the  basis  of  the  para¬ 
meters  of  e1  and  p=p  ,  the  parameters  of  point  c  can  then  be 

c 

calculated  from  formulas  (4)  and  (5) .  As  a  result,  the  nozzle 
exit  temperature  and  molecular  weight  formulas  are  obtained 


T, 


>  -  (rn 
M.  -  (m„ 


0j4n  -te-'j  exp  ^  Dr„  +  £Wa-^j|lii  - In  *  ^(12) 

5o»*ln-^-Vln  —  la  ^1  (13) 

friA  ?«  ,  *~/J  .  ., 


+  w.  £) ««.  [(. 


''Mm 


3.  Formula  of  Specific  Thrust  Parameters 


By  substituting  H  =H  ,=H  +  AH  and  H  “Hc^H  +Q  In 

C  C  CO  C  C  CO  H  p 


c1 


in  the  design  altitude  specific  thrust  formula 
1=294.98  -y(Hc-Hc)/1000,  we  obtain 


co 


/  -  9.328  Jfit.  -+•  AH,  —  —  QJa 

1  ?— 


9.:0  O 


The  interpolation  coefficients  QT#  Q, 


DT'  °M  'V  Md  °H  in 


formulas  (12)  to  (14)  use  values  in  the  vicinity  of  p  .  When 
_  _co 

Pc •  >  Pc • 

In —  <0,  values  smaller  than  p  are  used;  when  ln—^— >0  values 
d  co  p  ' 


co 


co 


are  greater  than  pcQ  are  used. 


The  nozzle  exit  ratio  area  and  nozzle  area  ratio  are  calcul 
ated  by  the  following  formulas; 


U  -  86.50 T,/lp,M, 

8  T.±*tfJc * 


05) 

(16) 


V.  ACCURACY  OF  EXTRAPOLATION  FORMULAS 


Extrapolation  calculations  can  be  carried  out  according  to 
the  above  formulas.  The  results  of  a  number  of  extrapolation 
calculations  for  oxygen-hydrogen  propellant  with  a  mixing  ratio 
0/F*7.9365,  the  nozzle  area  ratio,  characteristic  velocity, 
combustion  chamber  temperature  and  design  altitude  specific 
thrust  are  given  n  Tab)'  2,  3  and  Fig.  2  respectively.  For 
comparison,  the  cor'  ^sponding  results  of  exact  calculations 
[5,6]  and  the  results  of  extrapolation  calculations  from  refer¬ 
ence  5  are  also  given. 


From  tables  2-4  and  fig.  2  it  can  be  seen  that  the  accuracy 
of  the  extrapolation  formulas  of  this  paper  is  relatively  high 
and  can  be  suitably  used  for  extrapolation  calculations  of  the 
actual  parametric  range  of  rocket  engines. 

Table  1.  Parameters  on  isentropic  reference  line  (gas  hydrogen 
gas  oxygen  propellant,  mixing  ratio  0/F=7.9365  initial 
enthalpy  H  =0  Kcal/Kg) . 


W.  *+/4M f  (6) 


(7) 


u.  *-*/**«  (8) 


rm 

™( 

h _ 

K 

1.727 

m ■ 

U.I5 

39.40 

(.MI 

MM 

E 

Mil 

0 

-247.7 

-942.0 

n.m 

11. 774 

l(.(M 

2.54(1 

2.2704 

<0.040 

D 

n 

El 

□ 

El 

n 
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0.2264 

0.1701  i 

(.1114 

0.0091 

2409 

2210 

2131 

1909 

'HHIWWPyWBMft.Biti  wm&\  ftJWTTffgTrTBTffWHyi? 


17.11# 

17.421 

17.702 

17.7(1 

17.04 

17.(71 

17.707 

IZ3 

1.75M 

1.(277 

1.2(74 

1.17(7 

1.0177 

0.77(4 

0.7217 

ES& 

0.014(4 

0.01214 

0.00(11 

O.OOSQO 

0.00130 

0.M270 

» 

0.M212 

E3 

Key:  (1)  Parameters;  (2)  Combustion  chamber;  (3)  Throat; 
(4)  Nozzle  exit;  (5)  kg/cm2  (absolute) ;  (6)  Kcal/kg; 
(7)  kg/kg  mole;  (8)  Kcal/kgK;  (9)  Data  obtained  by 
extrapolation. 


(3).** 


Fig.  1  Enthalpy-entropy  schematic 
diagram. 

Key:  (1)  Enthalpy; 

(2)  Isentropic  refer¬ 
ence  line;  (3)  Entropy. 
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******  (5) 


»  M  MM  100  4M  1000 

(7) 


Fig.  2  Comparison  of  design  altitude 
specific  thrust  calculated  by 
extrapolation  with  exact  value 
(hydrogen-oxygen  propellant,  mixing 
ratio  O/F-7.9365,  combustion  chamber 
pressure  p  =6.805  kg/cm^ , 
absolute) . c 

Key:  (1)  Design  altitude  specific 
thrust  I  seconds;  (2)  Gas 
hydrogen  ozone  (enthalpy  H  - 
629.1  Kcal/kg) ?  (3)  Liquid0 
hydrogen  liquid  oxygen  (enthalpy 
Hc=190.6  Kcal/kg);  (4)  Exact 
value;  (5)  Extrapolation  value 
of  this  paper;  (6)  Extrapolation 
value  of  reference  5; 

(7)  Pressure  ratio. 


Comparison  of  nozzle  area  ratio  calculated  by 
extrapolation  with  exact  values  (hydrogen- 
oxygen  propellant,  mixing  ratio  0/F*7.9365 
combustion  chamber  pressure  pc=6.805  kg/cm2 
absolute) 


Iff*#*,  »<t  *«  -  -190-6  *+/& J?  -  629.1  X-t/ZfT 


3)11*41 


2.468 

7.151 

10.75 

14.69 

36.28 

46.04 

64.44 

81.77 

98.31 


Key:  (1)  Liquid  oxygen  liquid  hydrogen,  enthalpy  Hc=190.6 
Kcal/kg;  (2)  Gas  hydrogen  ozone,  enthralpy  Hc=629.1 
Kcal/kg;  (3)  Exact  value;  (4)  Extrapolation  value  of 
this  paper;  (5)  Extrapolation  value  of  ref.  5; 

(6)  Exact  value;  (7)  Extrapolation  value  of  this  paper 
(8)  Extrapolation  value  of  ref.  5. 


{1';.  *****  * 


)***«,**  » 

/  _ tOA  l  I  4 


9)* 


af*  1 ** 5 


»  *  *(14) 


■X  X*5 


2156  2156  2157  2324  .  2326 


2142  2144 


2136  2132  2139  2099  2094  2099  2251  1  2254 


(17)  ft: 


Table  3.  Comparison  of  characteristic  velocity  calculated 
by  extrapolation  with  exact  values  (hydrogen- 
hydrogen  propellant,  mixing  ratio  O/F-7.9365). 

Key:  (1)  Combustion  .chamber  pressure  pc,  kg/cm2,  absolute; 

(2)  Gas  hydrogen  gas  oxygen  enthalpy  Hc=0  Kcal/kg; 

(3)  Exact  value;  (4)  Extrapolation  value;  (5)  This 
paper;  (6)  Reference  5;  (7)  Liquid  hydrogen  liquid 
oxygen,  enthalpy  ^=-190.6  Kcal/kg;  (8)  Exact  value; 

(9)  Extrapolation  value?  (10)  This  paper;  (11)  Refer¬ 
ence  5;  (12)  Gas  hydrogen  gas  oxygen,  enthalpy 
Hc=629.1  Kcal/kg;  (13)  Exact  value;  (14)  Extrapolation 
value;  (15)  This  paper;  (16)  Reference  5;  (17)  Note;  in 
this  table,  the  unit  for  the  characteristic  velocity 
values  is  m/sec. 


Table  4.  Comparison  of  combustion  chamber  temperature 
calculated  by  extrapolation  with  exact  values 
(liquid  hydrogen  liquid  oxygen  propellant, 
mixing  ratio  O/F-7.9365) 


HRftSft  fc* 

(1) 

1.02 

2.04 

4.0* 

10.21 

20.42 

40.83 

61.23 

81.56 

(2) 

3039 

3127 

3217 

3341 

3437 

3534 

3591 

3632 

(3) 

3033 

3138 

3223 

3344 

3437 

3333 

3390 

3631 

(J_)  ft: 

2 

Key:  (1)  Combustion  chamber  pressure  p  kg/cm  , 

absolute;  (2)  Exact  value;  (3)  Extrapolation 
value  of  this  paper;  (4)  Note:  in  this  table, 
the  unit  for  the  combustion  chamber  temperature 
values  is  K. 
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A  THERMODYNAMIC  CALCULATION  METHOD  OF  THE  SYMBOLIC  FORMULA 
FOR  ROCKET  ENGINES 


Fang  Zhaokui 
Abstract 

A  symbolic  element  is  substituted  for  chemical 
elements  which  cause  practical  chemical  reaction.  With 
the  use  of  the  symbolic  formula,  a  new  calcualtion  method 
has  been  set  up  to  calculate  theoretical  performance  of 
the  propellant.  It  will  be  able  to  devise  a  general  com¬ 
puter  program  for  chemical  equilibrium  of  rocket 
engine,  which  does  not  depend  on  practical  chemical 
element. 

In  the  past,  a  chemical  thermodynamic  calculation  method  was 
widely  used  in  China.  Its  special  feature  is  that  the  chemical 
reaction  formula  must  be  correlated  with  the  practical  chemical 
elements  of  the  propellant.  It  is  therefore  very  difficult  to 
design  a  general  program.  The  symbolic  thermodynamic  calcula¬ 
tion  method  proposed  in  this  paper  eliminates  the  above 
difficulty. 

I.  SYMBOLIC  ELEMENTS  AND  THEIR  MOLECULAR  FORMULAS 


We  let  a  pair  of  positive  integers  i  and  j  produce  i  types 
of  propellants  from  combustion  products  for  any  element  having 
j  to  define  a  group  of  symbolic  elements  as  A^  and  their 
corresponding  symbolic  element  atomic  number  by  a ^ .  The  ordered 
permutation  of  symbolic  element  A^  is: 


JI  Apa  “  (»  “  1 ,  2 , ■  •  •») 


This  expresses  the  molecula  formula  for  i  types  of  combust.-’ 
products  produced  after  propellant  combustion.  The  symbol  ,  jj 
expresses  the  ordered  permutation  as  in  formula  (1) . 


For  example,  for  the  four  element  propellants  C,  H,  0  and 
N,  A^,  A2 ,  A^ ,  A^  respectively  are  used  as  symbolic  elements  in 
their  place.  Letting  i  types  of  combustion  products  be  CO, 
obviously,  on  the  basis  of  this  definition,  its  molecular  form¬ 
ula  for  the  symbolic  elements  is: 

i 

CO  -  n  "  AllAiQAllA,Q  _  (2J_ 

The  relationship  between  the  designational  order  i,  j,  and  the 
products  and  elements  is  arbitrary.  However,  as  soon  as  a 
relationship  is  assigned,  formula  (2)  only  determines  CO. 


II.  A  CHEMICAL  EQUILIBRIUM  EQUATION  FOR  SYMBOLIC  ELEMENT 
COMBUSTION  PRODUCTS 

1 

Definition:  For  any  molecular  combustion  product  U  , 

let  the  dissociation  reaction  of  the  correlation  of  the  product 
and  its  symbolic  element  A^ 

n  x  «*<<<  O'  -  u  2-  •  *0  o) 

be  a  formula  of  a  chemical  equilibrium  reaction.  Because  several 
partial  factors  are  assigned  to  the  symbolic  elements,  their 
corresponding  conventional  rules  of  multiplication: 


X  A, 


If  au  -  0 
if  ,«**o 


(4) 


By  the  theory  of  chemical  equilibraium  reaction,  the  cor¬ 
responding  chemical  equilibrium  equation  for  formula  (3)  will  be: 

1  X  PX  X  •••  XP^)-/,(T)  0  -1,2...,)  (5) 

P  P 

In  the  formula,  .  nv<<'  **'  respectively  express  the  partial 
pressure  of  the  symbolic  element  combustion  product  n  Aj*H 
and  the  atomic  product  A^  during  the  combustion  process,  f i (T) 
is  the  equilibrium  constant  of  the  dissociation  reaction  under 
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temperature  T.  If  a  certain  element  does  not  exist  in  this 
molecule  (e.g.  there  is  no  hydrogen  in  the  CO  molecule) ,  then 
because  of  the  conventional  mathematical  law  *1,  formula 

Aj 

(5)  will  inevitably  be  satisfied.  From  the  left  and  right 
sides  of  formula  (5) ,  the  following  logarithm  can  be  obtained: 


In  P, 


sl  ~  * 


By  using  formula  (6)  ,  molecular  partial  pressure  can  be 
converted  into  atomic  partial  pressure,  molecular  variables  can 
be  eliminated  and,  furthermore,  we  can  greatly  decrease  the 
number  of  equations. 


III.  PROPELLANT  MATRICES  FOR  SYMBOLIC  ELEMENTS  AND  THEIR 
NORMALIZATION  PROCESS 


Definition:  If  there  is  a  type  of  complex  propellant 
composed  from  n  types  of  pure  chemical  mixtures,  each  type  is  a 
pure  chemical  molecular  expression.  When  the  percentage  of  heat 
content  and  mixture  weight  is  respectively  JJ  x.  (-*  —  i,  2„  . 

then  the  characteristics  of  the  propellant  are  completely 
determined  by  the  matrix 


(7) 


Formula  (7)  is  called  the  propellant  matrix.  In  performance 
computations,  it  is  only  necessary  to  consider  the  theoretical 
performance  possessed  by  the  unit  mass. 

For  this  reason 


L4 


Definition:  The  unit  mole  propellant  symbolic  element 
atomic  formula  which  takes  the  molecular  weight  as  1  is  the 
standardized  equivalent  formula  of  the  propellant.  The  process 
of  converting  the  original  propellant  into  an  equivalent  formula 
for  a  symbolic  element  propellant  is  known  as  the  symbolic 
element  propellant  normalization  process. 

Theorem:  If  there  is  a  complex  propellant  (7)  ,  ju.^  is  the 
atomic  weight  of  the  symbolic  element  A^ ,  letting  the  matrix 


then  the  normalization  process  of  propellant  (7)  is  realized  as 
the  matrix  operation  of  the  following  formula 


*  »  '*  ‘  '  A  '  V  -  '  *  ' 


„  t  -  -  f  it  '  1 


Mj+O  *  >  ■*  ?(»)  >  3 

•  r.« :v.rm 

■  Z  1  '/f*'  •*».  n:7~:- £**£*«' 


.  f,  (9) 


■1 


In  the  formula.  A' (n, j+1)  is  the  transposed  matrix  of  A(n,j+1) 
obtained  after  taking  out  the  final  column  from  matrix  (7) . 

(9)  is  a  column  matrix  after  normalization.  The  subscript  "0" 
of  the  column  matrix  element  a^  shows  that  it  does  not  belong 
inside  the  combustion  product  serial  number  "i".  From  (9)  we 
can  deduce: 


‘i  :  hobmtO  r  ”*.1.* 


•  *. 


do) 


This  is  the  inevitable  result  of  normalization  processing.  For 
the  double  element  liquid  propellant,  if  the  incendiary  agents 
corresponding  to  the  forms  in  (7)  and  (8)  are  A(n,j+2) ,  B(n), 
the  oxidant  is  A(m, j+2) ,  B(m) ,  and  K  is  the  mixing  ratio. 
Therefore,  the  normalization  process  of  this  type  of  complex 
double  element  propellant  is: 

'  A(j  +  1)-  *  B(j *)  +  KA'(m,  i  .+  i>*  *C*XJ/0  +  K> 

The  final  result  is  similarly  a  column  matrix.  (9) ,  (10)  and 
(11)  are  proof. 

IV.  A  BASIC  EQUATION  AND  PROGRAM  DESIGN  FOR  THE  THERMODYNAMIC 
CALCULATION  OF  THE  SYMBOLIC  FORMULA 

Given  the  propellants  determined  by  the  positive  integers 
i  and  j  and  the  normalized  matrix  A (j+1) ,  the  standard  formula 
for  the  combustion  reaction  is: 

i  i  i  _  ...*  *  —• 

m  xi  E  .(,z> 

3  5t- 

The  above  formula  expresses  that  the  combustion  product 
has  partial  pressure  separately  produced  by  M  mole  propellants 
at  the  time  of  combustion.  However,  when  there  is  a  condensed 
phase,  P^  should  be  seen  as  the  number  of  moles  of  the  condensed 
phase  product.  The  four  basic  equations  for  the  symbolic 
thermodynamic  calculation  are  now  changed  to  the  following: 

1.  Equation  of  Mass  Conservation: 

I,2,***/  (13) 

2.  Equation  of  Energy  Conservation: 


H9  -  2  WXPi/M  (1) 

(**MMf#£)  (2) 


Key:  (1)  Isenthalpy-combustion 
chamber  state;  (2)  isentropy- 
nozzle  state. 


3.  Equation  of  Chemical  Equilibrium: 

,  (F*V(I?T)  4*ln-»-2  +1“^] -0  (15) 

In  formulas  (14)  and  (15),  (H®)^  (S°)  i  and  (F°)  are  respec¬ 
tively  the  enthalpy,  entropy  and  free  energy  of  i  type  products 
at  temperature  T. 


R\aPi  ^for  9as  phase) 
o  (for  condensed  phase) 


-fi 


pt  (for  gas  phase) 

l  (for  condensed 
phase) 


R  is  a  mole  universal  gas  constant.  The  discriminants  of  the 
solidifying  phase  are  omitted. 


4.  Dalton's  Law: 


(16) 

'  *  -t.r  ». 


Key:  (1)  Only  for  gas  phase, 


When  A(n,j+2),  A(m,j+2),  K,  PQ,  HQ  and  are  already 

known,  a  basic  proposition  of  a  thermodynamic  calculation  is 
formed  in  seeking  a  solution  to  the  above  equation. 


In  program  design,  the  system  of  equations  must  be 
linearized,  and  when  using  the  reduced  number  of  equation  (6) 
the  partial  derivative  approach  and  divergence  factors  are 
used  to  solve  the  problem  of  initial  value  arbitrariness 
within  a  wide  range.  We  used  BCY  language  to  work  out  a 
general  program  for  thermodynamic  calculations  of  the  symbolic 
formula  and  were  successful  in  calculating  logarithms  on  a 
109  C  machine  for  several  tens  of  types  of  propellants.  The 
above  contents  are  limited  in  space  but  further  details  are 
not  necessary. 
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A  CORRECTED  CALCULATION  METHOD  FOR  ISENTROPICALLY  EXPANDING 
PROCESSES  OF  COMBUSTION  PRODUCTS  OF  HYDROCARBON  WITH  AIR 


Yang  Donghua 

East  China  Institute  of  Chemical  Technology 


Abstract 


In  this  paper  an  immediate  and  more  accurate  calculation 
method  for  isentropically  expanding  processes  of  combustion 
products  of  hydrocarbon  with  air  is  suggested.  The  accuracy 
of  this  method  suggested  by  this  paper  is  4-9  times  more 
accurate  than  the  /Sr  number  method  [1]  and  the  B 

w  ri  2 

number  method  [1]  while  its  computation  is  a  little  more 
laborious. 


Three  types  of  approximate  calculations  for  isentropically 
expanding  processes  of  non-standard  combustion  products  are 
presented  in  reference  [1] .  A  simpler  and  more  exact  calculation 
method  for  isentropically  expanding  processes  of  non-standard 
combustion  products  (limited  to  hydrocarbon  combustion  products) 
is  proposed  in  this  paper. 


If 


-  L  ■ ••  fl-  /f*  •>  i,  •  f  lit 

and 


are  introduced,  then  the  non-standard  combustion  product  is 
isentropically  expanded  from  initial  state  1  to  final  state  2. 

Its  heat  drops  Al  and  Alqir0  can  be  respectively  corrected  as: 

n/  -■  J  CfJT  ™  <p(a/)ch,  Si 

HA725^T_^a  ..) 

J»  T 


Key :  ( 1 )  And . 


After  coefficients  *P  and  &  have  been  corrected,  it  is 
convenient  to  directly  use  reference  2  and  borrow  the  following 
steps  to  carry  out  the  calculation: 


Alg*'  — -  -  (Alg  *")/(/>  — * 

,  _ (*?/»J)ch,  —  consult  ref.  2  for  calculations^ 

Tit  (a/)ch,  ' — ►  A /  —  <p(A/)ch,  •  However ,  if  the  above  steps  are 

used  directly  from  reference  2  to  carry  out  the  calculations,  it 
is  first  necessary  to  have  calculated  the  corrected  coefficients 
and  for  the  gases  of  the  various  combustion  products  and 
numbers.  However,  for  reasons  given  below,  it  is  not 
necessary  to  first  calculate  the  corrected  coefficients  for  the 
gases  of  the  various  combustion  products  and  0  numbers  one  by 
one.  This  is  because  all  the  gases  of  the  hydrocarbon  combustion 
product  CHn  with  various  /9  numbers  can  be  viewed  as  being 
formed  from  the  gas  superposition  of  the  C  and  CH_  of  the  same 
fi  number: 

(1  +  0.25*  +  0.052476*/?)  (gas)  *ch,  “  O'—  0.5»)  (gas)  +  (0.75* 

■+■  0.052476*0),  ( gas )  ftCH, 

Therefore,  we  can  obtain: 

(gas)  fen*""  lA'/C.*'  +  B’’)].  (gas)  *c  +  [B'ftA'  ■+•  SOJi  (gas)  #»ch. 

In  the  formula,  A'+l-0.5n;  B ' =0 . 75n+0 . 052476n/ff  .  In  this  way, 
after  undergoing  isentropic  expansion,  the  corrected  values 
Cp  ^  and  H  of  the  gases  of  CHn  can  be  calculated  by  using  the 

corrected  values  and  of  the  gases  of  C: 

<pi  -  (A/f/n/J)  -  U'tf  +  BO/ U  +  BO, 

*1  -  C(Alg*S')/(Alg*P)]  -  ( A'rt  +  BOKA'  +  BO 


Thus,  provided  we  have  the  corrected  values  and  jJ ^  of  all 

the  gases  of  the  numbers  of  C,  the  isentropic  processes  of 
CHn  hydrocarbon  with  air  combustion  products  can  be  directly 
calculated  according  to  the  steps  described  above  and  using  the 
most  exhaustive  and  up-to-date  table  of  the  thermodynamic 
properties  of  gases  [2] . 

The  author  calculated  the  /£?  numbers  as  0.1,  0.2 . 1.0 

of  the  gases  of  C.  Ten  types  of  initial  temperatures  from 

1500,  1400 . 600 °C  etc.  were  isentropically  expanded  to  the 

corrected  coefficients  and  of  the  ten  final  tempera¬ 
tures  1050,  850 . etc.  From  the  results  of  the  calculations, 

it  can  be  seen  that  changes  are  slow  for  P^and  values  in 

relation  to  initial  and  final  temperatures.  Therefore,  this 
type  of  detailed  table  of  corrected  coefficients  can  be  used 
satisfactorily  for  engineering  calculations  [3]. 

We  later  carried  out  an  accuracy  analysis  of  the  calcula¬ 
tions.  From  Table  1  it  can  be  seen  that  the  accuracy  of  the 
corrected  calculations  in  this  paper  is:  heat  drop  Al  $  0.05%; 
final  absolute  temperature  T2  $  0.1%.  The  accuracy  was  5-10 

times  greater  than  the  fir  number  method  and  the  BCH  number 

2 

method,  although  the  calculations  were  slightly  more  complex. 
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Table  1  Using  the  corrected  calculation  method  in  this  paper, 
the  following  were  calculated  from  each  kind  of  initial 
temperature  t^°C:  (gas)  ,^=1, H2 ,  (gas)  ^=1  ,CH4 ,  (gas)/5>=0 . 2 , 

CH3*887'  /^=0-2,H2,  as  well  as  the  t2°C  and  values 

and  their  relative  errors  <^I%,  T2%  after  isentropic 

expansion  of  the  fixed  pressure  ratio  p^/p2 • 

Key:  (a)  Combustion  product;  (b)  True;  (c)  Content;  (d)  Pressure 
ratio;  (e)  True  Kcal/kg-mole;  (f)  True  Kcal/kg-mole; 

(g)  True;  (h)  True;  (i)  True;  (j)  True;  (k)  True. 
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THE  EXPERIMENTAL  INVESTIGATION  OF  THE  TIP  FLOWS  IN  DIFFUSER 
CASCADES 


Lai  Shengkai,  Su  Rongpei,  Su  Jiexian,  Cheng  Le  and  Jiao  Deyong 
Harbin  Institute  of  Technology 

Abstract 


The  results  of  the  experimental  investigation  of  tip 
flows  in  diffuser  cascades  are  presented.  The  results  of 
measurements  of  secondary  flow  velocities,  main  flow 
velocity,  total  pressure  as  well  as  air  flow  deflection 
angle  in  the  vicinity  of  the  end  wall  of  the  exit  plane  of 
the  cascade  are  presented.  The  authors  point  out  that  the 
effects  of  the  shrouding  on  the  aerodynamic  characteristics 
must  be  considered  when  inner  shrouding  is  employed  in 
compressor  stationary  cascades. 


In  order  to  determine  the  effect  of  movement  of  the  end 
wall  of  an  axial  flow  compressor  stationary  cascade  on  secondary 
flow,  the  authors  carried  out  experimental  investigation  on  a 
stationary  cascade  wind  tunnel  (M  no.  approx.  0.2).  The 
experimental  cascade  consisted  of  circular  arc  straight  vanes. 

In  order  to  simulate  an  unshrouded  cascade,  a  flat  variable 
speed  moving  belt  was  substituted  for  the  end  wall  surface  of 
the  experimental  cascade.  A  five-hole  probe  was  used  to  measure 
air  flow. 


Figures  1  and  2  each  show  the  distribution  of  secondary 
flow  velocity.  In  the  figures,  W  is  the  component  velocity 
along  the  direction  of  the  cascade  spacing;  the  direction  toward 
the  suction  surface  is  positive;  is  the  component  velocity 

along  the  direction  of  the  blade  height,  and  because  of  the  end 
wall  the  direction  toward  the  center  is  positive.  Wx  and 
in  the  figures  are  both  relative  values  of  the  main  flow 
velocity  W2q  for  the  same  cascade  spacing  locations,  v  is  the 
velocity  of  movement  of  the  belt. 


From  the  velocity  distribution  in  figs.  1  and  2  it  can  be 
seen  that  the  secondary  flow  is  more  or  less  a  triangular  shaped 
circulation  region  in  the  air  passage.  The  effect  of  end  wall 
movement  on  circulation  is  mainly  produced  within  a  relatively 
small  region  near  the  wall  corner  of  the  suction  surface  and  the 
flow  effects  are  small  at  some  distance  from  the  wall  corner. 

When  the  end  wall  moves,  the  flow  in  the  positive  direction 
near  the  end  wall  is  markedly  weakened  and  the  circulation 
region  is  reduced. 

Figure  3  shows  the  distribution  of  mean  flow  direction 
vortices.  Counterclockwise  circulation  is  positive.  The  symbols 
"  #  "  and  "  in  the  figures  represent  the  vortex  strength  when 

the  lower  end  wall  of  the  cascade  is  stationary  and  when  moving 
respectively.  It  can  be  seen  from  the  figure  that  when  the  end 
wall  is  not  moving,  the  vortex  is  strongest  at  area  1  near  the 
wall  corner  of  the  suction  surface.  When  it  is  moving,  there  is 
a  general  drop  in  vortex  strength  and  the  strongest  vortices 
shift  from  area  1  to  area  2. 

From  fig.  4  it  can  be  seen  that,  apart  from  near  the  suction 
surface,  the  effect  of  the  end  wall  movement  on  main  flow  velocity 
is  very  small.  Near  the  wall  corner  of  the  suction  surface  there 
is  a  marked  drop  in  main  flow  velocity  when  the  end  wall  is  not 
in  motion.  When  the  end  wall  is  in  motion,  there  is  a  relatively 
large  change  in  main  flow  velocity  near  the  wall  corner  of  the 
suction  surface  and  an  increase  of  main  flow  velocity  near  the 
end  wall. 

The  difference  between  the  air  flow  deflection  angle  of  the 
tip  region  and  that  of  the  main  flow  region  is  known  as  the  air 
flow  over-deflection  angle.  Figure  5  shows  the  distribution  of 
over-deflection  angles  (a)  and  the  mean  value  of  the  flow  rates 
(b) .  The  over-deflection  angle  is  greater  near  the  end  wall  of 
the  suction  surface  and  reaches  a  maximum  of  15°.  Because  of  the 


over-deflection  of  the  air  flow,  there  is  an  increase  in  blade 
loading  in  this  area  and  a  deviation  from  the  design  values  of 
cascade  inlet  angles.  Therefore,  the  characteristics  of  the 
tip  region  should  be  considered  in  calculating  the  design  of 
cascades.  End  wall  motion  can  cause  a  clear  reduction  in  the 
over-deflection  angle  in  the  wall  corner  region  of  the  suction 
surface. 

From  fig.  6  it  can  be  seen  that  flow  loss  is  mainly  con¬ 
centrated  in  the  area  of  two  wall  corners,  especially  the  area 
of  the  wall  corner  of  the  suction  surface.  End  wall  motion 
causes  a  reduction  in  high  loss  areas. 

Based  on  an  analysis  of  the  results  of  the  experiment  the 
following  conclusions  are  obtained: 

1.  There  is  a  lot  of  over-deflection  of  the  flow  in  the 
wall  corner  region  of  the  suction  surface  in  compressor  station¬ 
ary  cascades  because  of  the  presence  of  tip  flows. 

2.  End  wall  motion  can  weaken  secondary  flow  in  the  air 
passage  as  well  as  weaken  the  over-deflection  of  air  flow  in  the 
wall  corner  region. 

3.  To  know  whether  or  not  to  use  an  inner  shrouding  in  the 
stationary  cascade  in  the  design  of  compressors,  it  is  necessary 
to  make  adequate  allowance  for  the  effect  of  the  inner  shrouding 
on  air  flow  in  the  cascade. 


Fig.  1  Component  velocity  distribution  in 

in  the  direction  of  cascade  spacing. 

Key:  (a)  Blade  height  (%) ;  (b)  Cascade  spacing  (%) 


(a),i 


5 

U 

“™Se==S 

■ 

cT^ 

1 

IBM 

*  UJ  1 

hUs 

..  .. 

„ 

fTiSI 

*  1 

*— t 

iSdiisi 

(d) 


(X) 

- OJS 


•  0.2 


Fig.  2 


Component  velocity  distribution  in 
the  blade  height  direction. 

Key:  (a)  Blade  height  (%)  ; 

(c)  Pressure  surface; 


(b) 

(d) 


Cascade  spacing  (%) 
Suction  surface. 


Fig.  6  Total  pressure  isobars  of  tip  region. 

Key:  (a)  Blade  height  (%) ;  (b)  Cascade  spacing  (%)  ; 
(c)  Pressure  surface;  (d )  Suction  surface; 

(e)  Blade  height  (%) ;  (f)  Cascade  spacing  (%) ; 

(g)  Pressure  surface;  (h)  Suction  surface. 
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DATA  CORRELATION  FOR  THE  EFFECTIVE  FLOW  RANGE  OF  COMPRESSOR 
CASCADES 
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Abstract 


This  paper  presents  a  semi-analytical  method  for 
calculating  the  compressor  cascade  throat  width  and  the 
correlation  of  effective  inlet  flow  angle  range  with 
critical  mach  number  is  given.  Good  agreement  between 
the  prediction  and  test  results  is  observed. 

I.  INTRODUCTION 

Reference  1  summarizes  the  test  results  of  a  circular  arc 
mean  camber  line  C series  cascade  and  proposes  a  method  of  data 
correlation  of  the  inlet  flow  angle  which  changes  with  the  Mach 
number  and  is  determined  by  the  principle  of  twice  the  minimum 
loss.  A  suitable  design  for  the  angle  of  attack  is  chosen  by 
using  the  cascade  curve  a Mjjc  so  as  to  guarantee  the  high 
efficiency  of  the  design.  At  the  same  time  it  ensures  that 
there  is  adequate  stall  and  blocking  allowance  in  the  off-design 
operation.  These  design  requirements  can  be  satisfied  by  changing 
the  curve  ^  through  modifying  the  cross  section  geometry 
of  the  blade  row. 

When  modifying  the  compressor,  the  curve  a^^  M^  can  be 
used  to  analyze  the  operating  conditions  of  the  blade  row  cross 
section  so  as  to  determine  the  location  of  a  breakdown  or  of  a 
basic  level  of  inferior  performance.  A  program  of  rematching  is 
proposed  and  also  of  doing  quantitative  estimates.  It  is  there¬ 
fore  a  powerful  tool  in  guiding  compressor  design  and  modifica¬ 
tion.  It  is  of  immense  practical  value  to  apply  this  method  of 
correlation  to  our  commonly  used  A^g  compressor  profile. 
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II.  CORRELATION  METHOD 


The  correlation  method  is  basically  the  same  as  in  reference 
1#  except  some  of  the  correlation  curves  have  been  redrawn  on  the 
basis  of  data  in  reference  2  and  others  have  been  corrected  in 
accordance  with  the  difference  between  two  types  of  profiles 
(the  mean  camber  line  shape  and  the  equal  leading  edge  thick¬ 
ness)  . 


1.  Relative  Throat  Width  of  Cascade  o/s 

In  order  to  raise  accuracy,  the  method  in  reference  1  was 
not  used.  A  semi-analytical  method  of  calculation  has  been  used 
(see  fig.  1) . 

The  equation  for  a  parabolic  mean  camber  line  with  the 
maximum  curvature  at  45%  of  chord  length  is: 

.  ,  1/y  —  a g*,/*  +  <sg*j/(l  —  *); 

The  back  arc  coordinate  (x^,  y^)  can  be  found  by  the 
following  formula: 

.-'«5  —  7i  —~rl^^eoaai 

In  the  formula,  y^  is  the  half  thickness  of  the  profile 
and  ai  can  be  found  by  the  differential  of  the  mean  camber  line 
equation. 

m,  —  —ctgc^O  — 

In  this  way,  the  problem  of  finding  o/s  can  be  summed  up  as 

finding  the  minimum  distance  (o/c)  .  from  the  leading  edge 

_  min  3  J 

center  of  the  adjacent  blade  M(x  ,  y  )  to  the  back  arc. 

m  m 

Finally,  o/s*  [ (o/c) min“r/c] / ( s/c) ;  where  r  is  the  radius 
of  the  leading  edge  of  the  profile. 


2.  Maximum  Critical  Mach  Number  M^c 

The  inlet  angle  range  (or  angle  of  attack  range)  determined 
by  twice  the  minimum  loss  drops  with  an  increase  in  the  Mach 
number.  When  this  range  tends  to  be  zero,  the  relative  critical 
Mach  number  is  known  as  the  maximum  critical  Mach  number.  It 
mainly  depends  on  the  curvature  of  airflow  inside  the  blades  and 
on  blade  spacing.  By  using  the  data  in  reference  2  the  curve 
al/V/MNC  Can  1:36  ma^e'  ^nc  Can  k®  obtained  and  with  the  correlation 
of  the  relative  blade  angle  of  curvature  &  and  nodal  chord  ratio 
s/c,  the  curve  M^/J^  ^  2  **  ®  •s/c  can  be  obtained  as  in  fig.  2. 

In  the  figure,  are  t^le  correction  factors  taking 

into  consideration  the  effect  of  relative  blade  thickness 
t/c  (t/c=0.1  is  a  basic  reference  value)  and  the  inlet  blade 
angle  a^(a^=s45°  is  a  basic  reference  value)  on  M^.  The 
correction  curves  can  be  expressed  by  the  following  polynomials: 

-  1.03943  +  l.l0357(//c)  -  l4.6429(f/c)a 
<A,-  0.931  +  0.00157  •  a\ 

JJ 

3.  Optimal  Inlet  Angle  a^ 

The  corresponding  minimum  loss  inlet  angle  relative  to  the 

maximum  critical  Mach  number  is  called  the  optimal  inlet  angle. 

H 

For  a  profile  with  the  same  thickness,  a^  depends  on  o/s.  By 
using  the  data  in  reference  2,  it  is  also  possible  to  establish 
the  relational  curve  in  fig.  3. 

4.  Positive  and  Negative  Stall  Angle  of  Attack  Change  in 
Relation  to  Mach  Number 

The  relational  curve  between  the  change  in  the  positive 

H 

stall  allowance  (a. -a. ) /x1 *x,  and  negative  stall  allowance 


H  /"s  H 

(cos  a^-cos  a ^) / Q>  and  (M^-M^)  is  given  in  reference  1. 

Figure  4  uses  data  from  reference  2  to  obtain  two  relational 
curves. 

In  the  figure,  x^  takes  into  consideration  the  effect  of 

blade  thickness.  Because  the  relative  distance  between  the 

thickness  of  the  profiles  A4q  and  is  extremely  small,  the 

coefficient  x.  is  considered  to  be  identical  with  C„.  x_  takes 
1  4  2 

into  consideration  the  effect  of  the  nodal  chord  ratio  and  angle 
of  curvature  on  the  allowance  for  positive  stall.  By  using 
0=25°,  s/c=1.0  as  a  basic  reference  value,  the  curve  for  x2 
is  made  (fig.  5) .  &  takes  into  consideration  the  effect  of  s/c 

and  a^  on  the  allowance  for  negative  stall  (fig.  6) . 

III.  CALCULATION  RESULTS 

On  the  basis  of  the  correlation  method  introduced  in  this 

paper,  the  ai'v‘MNc  curves  for  14  cascades  from  reference  2  were 

calculated.  These  tallied  well  with  the  predicted  results. 

Using  numbers  4  and  5  cascades  as  examples:  the  difference 

between  the  a1  ^  curves  obtained  using  the  calculation  method 

in  reference  1  and  the  tests  was  relatively  large?  after  using 

the  method  in  this  paper  to  determine  o/s,  the  method  in  reference 

1  was  used.  Although  there  was  some  improvement  of  the  difference 

H  H 

between  the  experiment  and  the  a^,  the  positive  and  negative 

stall  boundaries  of  the  effective  flow  range  could  still  be 
observed  as  different  from  the  experiment.  The  third  type  clearly 
tallies  well  with  the  method  in  this  paper  (see  fig.  7) . 


VC*.,?.) 
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Fig.  1  Diagram  for  determining  relative 
throat  width. 


Fig.  2  Relational  curve  of  maximum 

critical  Mach  number  on  curvature 
and  nodal  chord  ratio. 


a - «•* 


Fig.  3  Relational  curve  of  optimal  inlet 
angle  on  relative  throat  width. 
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la)  --[i]«js 

(b)  -•  --iis  oil 

(c) 

(£)  a  Xttf2]4J»  No.5,  No.6 

Fig.  7  Comparison  of  calculation  results 
with  test. 

Key:  (a)  Method  of  reference  1;  (b)  Reference 

1  method  still  used  after  determining 
o/s;  (c)  Our  method;  (d)  Test  data  for 
nos.  5  and  6  cascades  in  reference  2. 
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THE  HODOGRAPH  METHOD  DESIGN  FOR  TRANSONIC  TURBINE  CASCADES 
ON  SURFACES  OF  REVOLUTION 

Chen  Zuoyi 
Qinghua  University 

Abstract 

Applying  the  basic  equation  of  hodograph  method 
corresponding  to  the  S.  surface,  the  calculated  methods 
and  equations  of  hodograph  design  for  transonic  turbine 
cascade  in  revolutionary  surface  are  presented,  and  on 
basis  of  the  simplified  condition,  the  solution  for  mixed 
type  equation  corresponding  to  the  compressible  function 
in  revolutionary  surface  is  derived,  the  calculated 
results  of  profiles  are  also  presented. 

Previous  hodograph  method  design  has  been  limited  to  plane 
flow.  This  paper  is  the  first  to  propose  velocity  coordinate 
flow  function  equations  corresponding  to  flow  surfaces  of 
revolution  and  integral  equations  returning  from  velocity 
surface  to  stream  surface.  These  serve  as  the  basic  equations 
for  the  calculation  of  the  subsonic  region  components  of  a 
transonic  turbine  cascade.  In  the  supersonic  region  we  started 
from  the  flow  function  equation  of  the  velocity  coordinate 
corresponding  to  the  surface  of  revolution,  and  derived  the 
compressibility  function  corresponding  to  the  flow  surface  of 
revolution.  Furthermore,  we  were  able  to  apply  mixed  type 
equations  to  the  surface  of  revolution. 

I.  FLOW  FUNCTION  EQUATIONS  OF  VELOCITY  COORDINATES  CORRESPONDING 
TO  THE  FLOW  SURFACE  OF  REVOLUTION  AND  THEIR  CALCULATION 
FORMULA 

Applying  the  basic  relationship  of  three  dimensional  flow 
to  velocity  potential  and  stream  function  and  the  basic  princi¬ 
ples  of  conversion  between  the  velocity  surface  and  the  physical 
surface,  we  deduced  the  velocity  coordinate  stream  function 


equation  corresponding  to  the  surface  of  revolution,  as  follows 


(1) 


In  the  equation,  b=b/bg  and  p  -  p/ p*  are  both  dimensionless 
quantities. 

For  the  cone's  surface  of  revolution,  ■VT  is  a  constant. 
By  taking  all  the  compound  derivatives  in  (1)  and  developing 
and  justifying  them,  and  then  inserting  the  regular  stream 
function  V^.,  L  (  ^r)  =L  (  ^) -L  (  )  is  satisfied. 

Therefore,  in  solving  the  equation  of  velocity  coordinate 
stream  function  equation  corresponding  to  the  cone,  it  is  only 
necessary  to  solve  the  partial  differential  equation  of  the 
following  regular  stream  function: 


.  The  relative  thickness  of  the  flow  blade  b  is  a  function  of 
M*,  9  .  Its  final  determination  depends  on  the  whole  hodograph 


design  of  the  three  dimensional  flow.  To  simplify  the  calcula¬ 
tions  in  this  paper,  we  used  the  linear  relationship  expression 
b«c^M*+c2  B  ,  where  c^  and  are  constants.  By  inserting  it  in 

equation  (4) ,  the  partial  differential  equation  of  the  regular 
stream  function  used  in  actual  calculation  is  obtained  as  follows 


In  order  to  fit  the  requirements  of  the  random  boundary 
node,  a  random  curve  grid  was  used.  If  M* , 0  is  an  orthogonal 
linear  grid,  m *}Q  is  a  random  curve  grid  and  the  included  angle 
of  M*  and  m*  is  a,  then  m*  and  B  are  equations  of  regular 
stream  functions  expressed  by  variables  as  in: 

>  ~  +'*-*** + + *  (*-•  db Sty _ 

•  ‘  A  -  ’•  *  -  • if.'  (fit- 7a  ■ 

~  — - f*  <1  —  fiilkr,  +  —  a  (4) 

'  coca  / 

*  .  .  ...  J  •  , 

In  the  equation,  the  coefficients  H^,  H2,  and  H4  are  all 
functions  of  M*  and  0  . 

II.  INTEGRAL  EQUATION  FOR  RETURNING  PHYSICAL  STREAM  SURFACE 

After  finding  J^r  by  equation  (4) ,  ^  can  be  found  on  the 
basis  of  the  relationship  between  ^ and  fce  .  The  profile 
coordinates  can  then  be  found  by  returning  the  physical  surface 
by  the  distribution  of  )/*  .  If  the  initial  radius  opposite  the 
cone  stream  surface  is  Rq  and  the  cone  angle  is  a^,  then: 

R  ™  +  Z 
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The  coordinates  of  the  surface  of  revolution  are  expressed 
by  Z,  *P  ,  and  R,  since  the  relationship  between  R  and  Z  on  the 
cone  has  already  been  determined.  For  this  reason,  in  finding 
the  profile  contour,  the  relationship  between  the  two  coordinates 
Z  and  is  all  that  is  required.  Their  integral  equations  are: 


,  ,  f*  kp-MZ  bo*8. 

*“*•  )*.  MV 


.  a(g)  d±  !  i  (ag) 

dd  dd  lpMm  d± 

dM* 

_  —  f*  hpmMZco*9m£nd 

V  <P*~k(R,  +  Ztg«t)M*p* 

.  8(g)a»  ,  i  (I wf 

dd  dd  ipM*  d<!> 

did* 


(d*V  ./1\ 

[de±  +  dB  (5) 

d<A  dM*  Ip  dM * 


Sills)  +  * 

dM*  Ip  dM* 


III.  THE  MIXED  TYPE  EQUATION  CORRESPONDING  TO  THE  SURFACE  OF 

VELOCITY  ON  THE  CONE  STREAM  SURFACE  AND  ITS  NOZZLE  SOLUTION 


Applying  a  nozzle  solution  to  the  supersonic  components 
during  design  means  there  is  little  change  in  the  contour  &  of 
the  turbine  cascade  and  relatively  large  change  in  M*.  If  we 
suppose  that  b=c^M*,  by  substituting  it  in  equation  (1)  and 
introducing  the  new  variable  dcr  =  -KdM* ,  we  get: 


—  M_  p&[—K(l>']  + 

«•  da  J 

P 

LM*  p  J 

Clearly,  if  equation  (7)  can  be  given  the  mixed  type 

equation  +  tygg  =0  of  the  plane  flow  by  making  the 

coefficient  of  %  zero,  then  to  satisfy  m*— —  — 

9  da  p 

K  must  be 


Equation  (7)  then  changes  to: 


(9) 


In  the  equation 

f  ]/<"•**>’ 

F(  C  )  is  the  compressibility  function  of  the  mixed  type 
equation  corresponding  to  the  surface  of  velocity  of  the  cone. 


If  the  expanded  Tricomi  approximate  solution  is  used  in  the 
design,  only  if  we  can  find  an  approximate  function  will  it  have 
the  form  of  an  expanded  Tricomi  compressibility  function  and  be 
near  enough  to  F(  &  )  in  numerical  value.  If  the  approximate 


function  F„  (  CT  )  -- 


s««r 


,  we  can  use  the  two  following 


a*  "  '  ;  ( 1-Cb  a ) 5 
boundary  conditions  to  determine  the  coefficients  B  and  C: 


,  \4Usi  1  wa  i 

L  da  J**-i  l  da 

2.  i  ■“  [F(v)]«*-i 

Furthermore,  the  coefficients  C  and  B  of  the  approximate 
compressibility  function  are  obtained. 

Figure  1  shows  the  profile  of  a  hodograph  method  design 
carried  out  on  a  cone.  The  design  parameters  are:  r=1.29, 

*0.238,  Q*,  =0,  T«*>  =1490K,  eout=730,  M*ut=1.15,  RQ=300mm, 

a^*26.56°,  HQ*50.05mm.  From  the  calculations  we  can  know  that 
the  difference  in  laws  of  velocity  distribution  and  similar 
parameters  is  extremely  large  on  cones  and  on  planes.  There¬ 
fore,  the  effect  of  three  dimensions  must  be  considered  in 
transonic  cascade  design. 


(10) 

(11) 
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Fig.  1  Hodograph  method  design  for  transonic 
turbine  profile  on  a  cone. 

Key:  (a)  Radian;  (b)  Millimeters. 

BASIC  AERODYNAMIC  EQUATIONS  EXPRESSED  BY  PARTIAL  DERIVATIVES 
ALONG  AN  ARBITRARY  STREAM  SURFACE  IN  SEMI -ORTHOGONAL 
CURVILINEAR  COORDINATE  SYSTEMS 


Wang  Zhongqi 

Harbin  Institute  of  Technology 


Abstract 


Basic  aerodynamic  equations  expressed  by  partial 
derivatives  along  an  arbitrary  stream  surface  in  semi- 
orthogonal  curvilinear  coordinate  system  have  been  derived 
in  this  paper.  These  equations  are  applicable  not  only  to 
relative  stream  surfaces  S1  (when  /^Tg^=l)  but  also  to  S2 

(when  g."  =r)  .  They  are  also  suitable  for  cylindrical 

1  ^  3  2  1 

coordinate  system,  if  x  =1,  x  =  -  and  x  =z. 


I.  BASIC  AEROTHERMODYNAMIC  EQUATIONS 


We  used  a  semi -orthogonal  curvilinear  coordinate  system  of 
corresponding  unit  base  vectors  and  by  using  formulas  (11)  -  (17) 
in  reference  1  we  were  able  to  express  the  vector  form  of  the 
basic  equation  of  aerodynamics  as  the  scalar  quantity  form  in  a 
curvilinear  coordinate  system. 

1.  Equation  of  continuity 

0x*  dx*  Qx*  ” 

2.  Equation  of  motion 

w  ,  w‘  rr 

Qx>  grr 

W  r  a  (<*/*>.)  9 (v^ y r)1  ,  WT  rr  ^ 

Qx'  a*-  J  Vi  Jin 

mm  9l  —  T  dt 
df  0X* 

In  the  equation,  the  corner  symbols  ,  jS 


(2) 

and  T 


are 


assigned  values  according  to  the  "cyclic  law."  Therefore, 
equation  (4)  separately  expresses  a  partial  equation  of  motion 
along  the  direction  of  three  curvilinear  coordinates. 

£2a(a=l,2,3)  expresses  the  inverse  variant  physical  component 
of  the  angular  velocity  vector  6)  of  turbo-machinery  revolutions. 
By  using  the  triangles  och,  oab  and  ofd  in  fig.  1,  it  is  not 
difficult  to  obtain: 

&  —  <a(cocfcsM  +  anfawfljjeofVS) 

OF—  —  mang/S 

&  ”  — «*(c ml  +  mag  cot  qcotdu/S)  *in  ^ 

V 

In  the  formulas,  s  *■  Vl  — 

3.  Equation  of  energy 

9l  +  ws  91  _ 

II.  BASIC  AERODYNAMIC  EQUATION  ALONG  A  STREAM  SURFACE 

The  relationship  between  the  partial  derivative  of  parameter 
q  along  the  stream  surface  and  the  ordinary  partial  derivative 
can  be  expressed  as 

Bq  9q  .  '-iVtu  dq  Bq  _  dq  NlSgn  dq  ,-x 

a**  a*1  **  N*/Tn  a** 

By  using  this  expression,  the  basic  aerodynamic  equation  expressed 
by  a  partial  derivative  along  the  stream  surface  is  obtained  from 
equations  (1) ,  (2)  and  (3)  respectively. 


(3) 

(4) 

(5) 


1.  Equation  of  continuity 


In  these  two  equations,  a  ,  /$  and  Y  are  assigned  values 
according  to  the  "cyclic  law.”  It  can  be  proved,  in  the  three 
partial  equations  of  motion  along  the  stream  surface,  that  only 
two  are  independent.  Therefore,  in  the  two  above  expressions, 
a  can  only  have  the  two  values  of  1  and  2  and  similary,  {3  and 
Y  can  only  have  two  values.  The  partial  equations  of  motion 


1  2 

along  the  stream  surface  in  the  x  and  x  directions  can  now  be 

obtained  from  equation  (11) .  Furthermore,  from  equation  (12) 

it  can  be  seen  that  f*  includes  two  components;  one  component  is 

a  3 

related  to  the  changes  in  force  in  the  x  direction  (changes  in 
pressure  and  centrifugal  inertial  force)  and  the  other  is  related 
to  changes  in  the  parameters  0^2»  ^ll'  g22  and  g33  in  the  x3 

direction  of  the  curvilinear  coordinate  system.  It  can  be 
proved  that,  if  f*(a=l,2,3)  is  taken  as  three  covariant  components 
of  a  vector  f*,  then  f*lw, 


3*  Equation  of  energy 

wx  3i  -  to1  di 


(13) 


It  can  be  seen  from  this  that  there  is  no  change  in  the  form  of 
the  equation  of  energy,  whether  in  coordinate  surface  x3=constant 
which  is  a  surface  of  revolution  or  in  two  classes  of  semi- 
orthogonal  curvilinear  coordinate  systems  of  a  meridian  surface. 


III.  SEVERAL  SPECIAL  CONDITIONS  OF  PARTIAL  EQUATIONS  OF  MOTION 


1.  When  coordinate  surface  x3=const  is  an  arbitrary  surface 
of  revolution  (  ^Tg^=l)  ,  equation  (11)  is  simplified  into  two 
partial  equations  of  motion  (a=l,2)  expressed  by  partial 
derivatives  along  the  stream  surface  S1 . 


w •  t£ 

_  5l_  _  _  5 1 _ 

3*-  3*-  U 


J 


+  +  2^/^Z  an du (W*ar—WT&') 

ox*  '  .  -  -  i 


(14) 


2.  When  coordinate  surface  x  =const  is  an  arbitrary  surface 

3 

of  revolution  and  is  also  a  stream  surface,  W  =W3=0,  f*=f*=0 

and  equation  (14)  is  simplified  into  equation  (47)  in  reference 
2.  It  can  be  proved  that  the  two  equations  given  in  (14) 
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(corresponding  to  a=l,2)  are  identical. 

3.  When  coordinate  surface  x^=const  is  a  cylindrical 

3  12 

surface  (x  »r*const)  and  x  =z,  x  =  - T  ,  we  obtain 

mm  90°,  Vfu  —  It  Vga~  —  Wt  —  PP'',,  WJ  —  PF,  —  —  P*^, 

"  wr*  Af,  “  *«»  AT,  “  —«♦»  Afj  “■  «„  —  0 

Under  these  conditions,  equation  (14)  is  simplified  into  equation 
(39)  in  reference  3. 

3 

4.  When  coordinate  surface  x  =const  is  a  meridian  surface,, 

we  obtain  Yg^-r ,  x  =  <p  ,  Q  *0 ,  W  =W3=W9  , rj  =0 ,  Z  |  =  Z 

and  by  using  formulas  (3)  and  (4)  we  can  obtain  equation  (92)  in 
reference  2  from  equation  (11) .  From  these  conditions 
dr/dy  <—  dcomd-o/dcp  —  d*fgvJd<p  —  d‘/g3Jd<p  —  0,  and  equation  (12)  changes 
into  equation  (93)  in  reference  2. 


Fig.  1  The  deduced  inverse  variant  physical 
component  a*  expression  of  angular 
velocity  to  is  used  to  chart  the  e^ 

and  e^  coplane,  and  the  CO  and  e^  coplane. 
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MEASUREMENT  OF  AERODYNAMIC  PARAMETERS  BEHIND  ROTORS" 


Liu  Sihong 

Northwestern  Poly technical  University 


Abstract 


This  paper  analyzes  and  discusses  the  measurement 
errors  which  occur  when  measuring  the  total  temperature 
and  total  pressure  behind  rotors.  Focus  is  placed  on 
solving  the  following  problems:  1.  The  error  equations  of 
rotation  state  are  derived  by  sweep  integration  for  the 
first  order  system  to  square  wave.  It  provides  a  relatively 
clear  theoretical  outline  for  the  measurement  errors  behind 
rotors  and  the  method  and  equations  to  calculate  the  errors 
of  rotation  state.  2.  The  mechanism  of  aerodynamic  blockage 
and  a  new  definition  of  aerodynamic  ratio  are  briefly 
presented.  An  experimental  curve  of  the  blockage  errors 
behind  rotors  and  the  ratio  of  aerodynamic  blockage  is 
given.  3.  The  method  to  decrease  the  measurement  errors 
behind  rotors  is  explored. 


I.  MEASUREMENT  ERRORS  BEHIND  ROTORS 

The  aerodynamic  parameters  of  total  temperature  and  total 
pressure  are  measured  behind  rotors  by  a  probe.  Because  of  the 
rotation  of  the  rotors,  it  is  possible  on  the  one  hand  to  have 
the  probe  measure  the  mean  operating  parameters  on  the  correspon¬ 
ding  flux  under  conditions  of  non-circular  motion.  On  the  other 
hand,  there  can  also  be  errors  in  measurement  on  account  of  the 
rotation.  These  types  of  errors  are  sometimes  positive,  some¬ 
times  negative,  sometimes  large  and  sometimes  small.  They 
cause  scattering  in  the  measured  values  and  result  in  distortion. 
We  analyze  the  tot' 1  temperature  as  an  example  below. 

1.  Rotation  State  Errors 


In  turbine  performance  experiments  it  is  the  mean  parameters 
of  the  ring  surfaces  and  the  whole  outlet  cross  section  which 
need  to  be  determined.  In  order  to  obtain  the  mean  total 

•This  paper  was  printed  at  the  Xiamen  Thermodynamic  Symposium 
in  August  1981.  -n 


temperature  of  the  rotor  outlet,  the  cascade  spacing  t_  must  be 

I, 

based  on  the  mass  mean 

TZ  -  j'L  T*C»pdz/\'  C„pdx  ( l ) 

If  it  is  incompressible  and  P  =const,  then 

TZ  -  \‘L  T'C^dx/^  C1Adx  (2) 

In  order  to  transform  the  space  field  into  a  time  field,  on  the 
basis  of  circumferential  velocity,  u=const,  x=ut,  dt=dx/u 
replacement  is  carried  out  on  equation  (2) ,  and  we  obtain 

T2  -  £ TTCuitf^C^t  >’■  (3) 

,  \  _ _ _ _  T  _ _  _  _ ... 

However,  the  thermocouple  behind  the  rotor  can  only  be 
affected  by  time  mean  values. 

TIT--fVw/,7-— '  (4) 

Because  the  mass  mean  requires  an  axial  velocity  C2A  weighted 

mean,  and  the  thermocouple  a  real  time  mean,  that  is  an  isobar ic 

mean,  along  with  the  T*  error  value,  it  is  an  error  introduced 

m 

with  unsuitable  weighting  called  the  weighted  error.  It  is 
expressed  by 

,  ^ (5) 

Key:  (1)  Weighted. 


When  the  rotor  sweeps  past  the  thermocouple  at  high  speeds, 
the  temperature  T ^  measured  at  the  thermocouple  contact  point 
cannot  keep  up  with  changes  in  the  total  temperature  of  tht' 
airflow  because  of  the  thermal  inertia  of  the  thermocouple. 
Therefore,  the  thermocouple  can  only  give  time  mean  values  of  T. 


»*»  s 


These  types  of  errors  introduced  because  of  inertia  are  known 
as  dynamic  response  erros  and  are  expressed  by  A T 

_ (1).  .  reSP°nSe 


Key:  (1)  Response. 

Weighted  errors  and  dynamic  response  errors  are  caused  by  the 

rotor  being  in  the  rotation  state.  Hence,  they  are  both  known 

as  rotation  state  errors  and  are  expressed  by  4T  .  ^ . 

v  2  rotation 

Mi  -  (2)  (2)  . 

-  iI-  +  AT.  „  .. gy ^  (») 

Key:  (!)  Rotation;  (2)  Weighted;  (3)  Response. 

There  can  also  be  these  types  of  errors  in  the  measurement 
of  total  pressure  behind  the  rotor. 


-  JV Cudt/^%CiAdt 


A  <■») 

A?.  ”  PH,  -PZ  -APm  +  &Pm  (11) 

(4)  (5) 

Key:  (1)  Weighted;  (2)  Response;  (3)  Rotation; 

(4)  Weighted;  (5)  Response. 


2.  Aerodynamic  Blockage  Errors 


The  probe  located  downstream  of  the  rotor  can  increase  the 
corresponding  back  pressure  of  the  channel  outlet  and  produce 
aerodynamic  blockage.  Aerodynamic  blockage  is  not  only  deter¬ 
mined  by  the  ratio  of  the  windward  area  of  the  probe  and  the 
channel  outlet  f/F^,  but  is  also  related  to  the  characteristic 
measurement  d  of  the  blockage  matter  and  its  distance  1  from  the 


channel  outlet.  For  this  reason,  if  the  aerodynamic  blockage 
is  only  measured  in  area  ratio  as  is  done  for  geometric  blockage, 
the  results  will  not  be  satisfactory.  Here,  fd/F^^  is  defined 
as  the  aerodynamic  blockage  ratio  and  it  is  used  to  indicate 
the  level  of  aerodynamic  blockage.  Figure  1  shows  the  experi¬ 
mental  curve  of  relative  blockage  errors  behind  the  double  stage 
turbine  of  a  turbojet  which  change  with  the  aerodynamic  blockage 
ratio.  The  results  shov,  that  the  blockage  causes  the  measured 
value  of  temperature  T*  behind  the  turbine  to  be  on  the  high 
side.  When  the  aerodynamic  blockage  ratio  reaches  16.5%,  it 
corresponds  to  the  blockage  error  of  the  maximum  state  M  and  the 
rated  state  H,  a  turbine  temperature  of  2.7%  and  2.9%  of  A TT 
respectively.  With  the  rearward  shift  of  the  probe  and  after 
the  aerodynamic  blockage  ratio  has  dropped  to  6%,  A  K. 

0.25%. 


Fig.  1. 

II.  THE  ANALYTICAL  ANALYSIS  OF  ROTATION  STATE  ERRORS 

The  flow  fields  behind  each  channel  of  the  rotor  are  not 
completely  identical,  however  an  equivalent  mean  characteristic 
can  always  be  used  as  a  substitute.  Consequently,  rotation  state 


errors  can  be  obtained  from  an  analysis  of  the  measured  values 
behind  this  channel.  In  view  of  the  fact  that  temperature  and 
pressure  signals  along  the  cascade  spacing  are  frequently  not 
continuous,  sweep  integration  is  used  here  to  derive  the  rotation 
state  error  equation.  So-called  sweep  integration  consists  of 
first  dividing  the  cascade  spacing  into  several  subintervals  and 
making  each  subinterval  satisfy  the  requirements  of  continuity. 
Afterwards,  there  is  successive  sweeping  beginning  at  the  first 
subinterval.  The  probe  is  then  placed  at  the  final  value  of  the 
preceding  subinterval ,  that  is,  it  enters  this  subinterval  as  an 
intitial  value,  and  there  is  one  by  one  integration. 

For  the  purpose  of  conciseness,  a  flow  field  behind  the 
rotor  approximating  a  rectangular  wave  fitting  was  used  (fig.  2) 
and  given  the  probe  is  a  first  order  system,  rotation  state 
errors  of  total  temperature  and  total  pressure  are  respectively 
discussed  under  conditions  where  the  trail  relative  thickness  b 
is  both  relatively  thick  and  relatively  thin. 

For  fig.  3a,  the  condition  is  shown  where  b  is  relatively 
thin  and  the  choice  of  the  coordinate  can  be  divided  into  the 
two  subintervals  0  ^  t^  and  t^^^.  Each  subinterval  is  a 
step.  The  step  quantity  for  the  total  temperature  is  D=T*-Tjg. 

Because  the  transfer  function  of  a  one  order  system  is 
G(s)*l/(1+Ts)  ,  therefore 

-fipi/cw 

The  thermocouple's  temperature  can  be  obtained  according  to  the 
Laplace  inverse  transform. 


We  substitute  formula  (12)  into  formula  (6)  and  integrate 
according  to  the  following  conditions 


a  «*£<«•  T*  -  Ti,  «  'v  *'•  •  » 

^<^,s  r*j*TTi  c^-r^r^  ?:,  +  (rr —t:)(i--  ,-**"). 


and  let 


3  ■»  tji,  “•<»/*/  *  */<£. 


*«- (Tt^TTXt-  '&*yt&*+'*i**y  - 1 1  (i# 

■-  - ,  -  Vi  ~- '  -  ..  ...  -  -  — <  ■*■  «  -  "*• 


Further,  from  formula  (13) ,  we  obtain 

-  •  -  \  •  ^  '  -  !«Jb  ;'*»«« ju  .  • 

Ti:**  (,TtSC»  _+  T«  3Cj^/ («<:*  -t^  iC)f^  ^  pj 


(i+J 


The  equation  of  rotation  state  errors  is  then  obtained  by 
substituting  formulas  (13)  and  (14)  into  formula  (8) 


AT,  -  3(1  -  S)(T*  —  T*)(C*  -  C,)/(3C*  +tcm) 

(1)  +  -  (77  —  1T?)(1  -  ,-‘»"X(*-'*"  _  ,-V)  -  U  j 


(15) 


Key:  (1)  Rotation- 


I. 


S 


For  the  same  reason,  when  b  is  relatively  thin,  the  equation  of 
the  rotation  state  error  of  total  pressure  is: 


Cl) 


—  V-W  -  |  •  :  t 

4£m.rJtLr  WS  -  F?)(c,.  -  <:»)/(«*  + 


^  - 


+  ^  (FT  -  F?)(l  -  _  e-'/»>  -■  i  J 

Key:  (1)  Rotation. 


(16) 


When  b  is  relatively  thick  (fig.  3b) ,  the  respective  equations 
for  the  rotation  state  errors  of  total  temperature  and  total 
pressure  are: 

_ m _  _  .  . .  ..  . . . 

AT,  -  3(1  -  5)(T»  -  T?)(C„  -  C1*)/(3C„  +  lC») 

+  —  (7?  -  T?)(l  -  e-'^'Xl  -  (*-*«"  -  e-S')]  (17) 

t, 

AF,  -  3(1  -  3 )(FT  -  F?)(C„  -  CJ*)/(5Cfc  +  JC,).  j 
(2) 

+  (pf  -  FT)(1  -  e-S")l  1  -  (*-'«"  -  «-**)  ]  (18) 

Key:  (1)  Rotation;  (2)  Rotation. 
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For  the  condition  shown  in  fig.  3,  we  know  from  equations  (5) 
and  (9)  that  the  respective  equations  for  weighted  errors  of 
total  temperature  and  total  pressure  are: 

(D  ...  -  -  -  . 

47*  -  i(l  -  5)(Tf  -  rJXC*  -  CJ*)/(5Clr+  *C*)  (19) 

4?."-  5(1  -  i)(Pf  -  P?)(C*  -  C1*)/(iCa.  +  iCu)  (20) 

(2) 

Key:  (1)  Weighted;  (2)  Weighted. 

Because  T*  >  T*,  P*  ^  P*/  a  ^  1  and  C2a  >  C2t),  it  is  not 

difficult  to  see  from  formulas  (19)  and  (20)  that  ^Twe!gjlte(j 

and  .  .  .  .  <*0.  This  is  mainly  because  the  trail  region 

weigntea 

flow  speed  is  small  and  the  weighting  corresponding  to  the  mass 

mean  time  is  relatively  small.  However,  the  mean  time,  i.e.  the 

isobar ic  mean  time,  is  actually  on  the  high  side  for  the  weighting 

of  this  section.  As  a  result,  under  the  condition  of  T*  >T*, 

b  ^  a 

T*  is  caused  to  be  on  the  high  side.  On  the  contrary,  under 
n 

conditions  of  P*  <  P*,  P*  tends  to  be  low. 

O  3i  n 

In  Korostelov's  report,  after  analyzing  a  triangular  profile, 
he  thought  that  when  weighting  was  unsuitable  the  opposite  con¬ 
ditions  would  appear,  i.e.  that  the  total  temperature  would  tend 
to  be  low  and  the  total  pressure  high.  This  clearly  differs 
from  the  results  discussed  above,  and  for  this  reason,  further 
analysis  was  carried  out. 

Korostelov,  in  analyzing  the  triangular  profile,  recommended 
the  following  conditions:  *  —  o~«,  *  —  o~«,  c-u— 

Ct,  P 0  —  Pt,  x  —  •  —  tt»  ft  Cii  ”  c*  —  °»  **  “  K  —  p*  +  Pc*/2* 

pCi/2,  P,  —  P*.  By  substituting  these  conditions  into 

formula  (20) ,  the  same  formula  as  in  the  above  paper  is  obtained. 

47,  -( 1  -  fXPt  —PT')  -  (1  —  «)(**  —  Cl  )p/2  '  V 

(1) 

Key:  (1)  Weighted. 

Because  a  (  1,  when  circular  velocity  is  u  >  C  #  it  can  happen 

<* 
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that  ^pweighted  ^•e’  the  total  pressure  value  tends  to 

be  high.  For  the  same  reason,  the  total  temperature  measured 
value  tends  to  be  low.  It  appears  that  this  conclusion  is 
correct.  Furthermore,  under  conditions  of  static  pressure 
equivalence,  Cfe=u  ^  Ca  which  implies  that  P*  >P*,  and  that  the 
total  pressure  of  the  trail  region  is  instead  higher  than  the 
main  stream  region's  total  pressure.  Naturally,  these  conditions 
could  actually  not  occur.  For  this  reason,  it  can  be  stated 
with  certainty  that  because  weighting  is  unsuitable,  the  total 
temperature  measured  value  behind  the  rotor  tends  to  be  high 
and  the  total  pressure  tends  to  be  low. 

Using  formulas  (19)  and  (20)  to  extract  the  partial  deriva¬ 
tive  of  b,  and  taking  into  consideration  a=l-b,  we  can  obtain 

€*)  +  (i  —  J)  KTf  -  TSXCa,  -  cu}  (21) 
dAP«/d5 -(?(<:,,--<:»)  +  («-*)](*? -tfX^-c*)  .  (22) 

(2) 

Key:  (1)  Weighted;  (2)  Weighted. 

It  is  thus  evident  that  the  weighted  errors  differ  with  b,  and 

-2  _  _ 

are  largest  when  b  (C2a-C2b)  =  (b-a)  . 

When  b  is  relatively  thin  (fig.  3a) ,  the  dynamic  response 
errors  of  the  total  temperature  and  total  pressure  are 
respectively: 

ATm  -  -  (7*  -  TTX 1  -  (<r-V*  _  ,-’✓»)  _  j  ]  (23) 

(1)  *'  r 

7 “  fc)(l  “  *-'»*)(  (*“'•"  -  r-'O  -  1  ]  (24) 

Key:  (1)  Response;  (2)  Response. 

When  b  is  relatively  thick  (fig.  3b) ,  the  dynamic  response 
errors  of  the  total  temperature  and  total  pressure  are 
respectively: 


5 


(25) 


CD, 

ATm  -  —  (Tf  -  r?)'(  1  -  -  ,-*✓*)] 

^  -  /?)(l  -  i  —  (,-»«"  -  <-/«)]  (26) 

12)  *'  - _  *  ;. 


Key:  (1)  Response;  (2)  Response. 


Because  «-v»  <  l <  ;;  i  >  (i  —  *-*V*) >  o,  l  >  (l  —  i-J*y  >  o;  -  [(«.-**"— 

‘  «-v<r)  —  1J  —  -  *-'**)  -  ll  <  o,  H  —  (*-v»  —,-*/»)]  -  {i  —  P-v»  (i  — 

«’“*/,)l  >  o;  .T? >  7?,  77  <Pf.  /  therefore  we  know  from  formulas  (23), 
(24) ,  (25)  and  (26)  that: 


When  b  is  relatively  thin, 
ture  tends  to  be  low;  when  A  P 
to  be  high. 


A  T 

response 

>0, 

response 


K.  0  and  the  tempera- 
the  pressure  tends 


When  b  is  relatively  thick,  4T  ^0  and  the 

r@spons6 

temperature  tends  to  be  high;  when  <0,  the  pressure 

irGsponso 

tends  to  be  low.  This  is  mainly  because  of  the  inertia  of  the 

« 

probe  which  can  neither  jump  nor  fall  sharply  with  it.  As  for 
the  total  temperature,  the  insufficient  part  of  the  jump  is 
greater  than  the  surplus  part  of  the  sharp  drop  when  b  is 
relatively  thin;  when  b  is  relatively  thick,  it  is  smaller  than 
the  insufficient  part  of  the  sharp  drop.  The  respective  values, 
therefore,  tend  to  be  lower  and  higher.  The  opposite  is  true 
for  total  pressure. 


The  formulas  also  show  that  the  dynamic  response  errors 
are  proportional  to  T/t  *f  /f  •  Therefore,  in  order  to  reduce 
dynamic  response  errors,  a  small  inertia  was  used;  the  probe  with 
the  relatively  large  cut-off  frequency  f*^  was  extremely  advanta¬ 
geous.  Because  t  *60/ (nz),  the  rotor  blade  number  increases  and 

3 

the  rotational  speed  accelerates.  As  (T*-T*)  and  (P*-P*)  enlarge, 
the  corresponding  dynamic  response  errors  also  increase. 

To  summarize,  in  carrying  out  measurements  behind  the 
rotor,  the  total  temperature  measured  value  tends  to  be  high  and 
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the  total  pressure  measured  value  tends  tends  to  be  low  because 
weighting  is  unsuitable.  Because  of  the  inertia  of  the  probe, 
total  temperature  is  lower  and  total  pressure  is  higher  when  b 
is  thin;  when  b  is  thick,  total  temperature  is  higher  and  total 
pressure  is  lower.  Rotation  state  errors  which  cause  the  meas¬ 
ured  values  to  be  high  or  low  should  be  concretely  analyzed 
according  to  actual  conditions. 


For  different  rotors,  if  their  (T*-T*)  and  (P*-P*)  are 

d  a  d  a 

relatively  large,  the  rotor  blade  number  is  relatively  large  and 
the  corresponding  errors  are.  relatively  large.  As  for  rotors 
which  are  the  same,  their  rotation  state  errors  will  not  only 
differ  with  the  position  of  the  radial  but  also  change  with  the 
rotational  speed  and  rate  of  flow. 

If  a  compressor  is  operating  at  high  efficiency  and  the 

probe  response  is  slow,  when  the  dynamic  response  errors  occupy 

a  main  position,  the  total  temperature  measured  value  will  be 

low  and  the  total  pressure  will  be  high.  However,  when  the 

compressor  is  close  to  stalling,  i.e.  when  b  is  relatively 

thick,  the  total  temperature  will  tend  to  be  high  and  the  total 

pressure  will  tend  to  be  low.  Figure  4  shows  two  methods  of 

measuring  the  rotor's  temperature  rise:  one  is  the  temperature 

rise  A determined  by  a  torsion  gauge.  It  can  approximately 

be  regarded  as  (T*-Tq) •  The  other  method  is  to  use  seven  ring 

surfaces  of  a  copper-constantan  thermocouple  parallel  to  the 

airflow  and  with  a  diameter  of  0.13mm  located  in  the  0.4b  area 

behind  the  rotor  to  determine  the  rotor's  temperature  rise  AT^. 

If  blockage  errors  are  neglected,  ^Tt*THj”Tg*  Their  error  value 

(4T  -4T  )  is  then  the  rotation  state  error  AT  .  .  . 
t  m  rotation 

Obviously,  ^Trotation  decreases  with  the  rate  of  flow  and 

increases  with  the  trail  relative  thickness  b.  The  laws  of 

change  for  the  gradual  increase  from  negative  to  positive  are  in 

complete  agreement  with  the  above  anaysis. 


0 


The  rotation  state  error  equation  cannot  only  be  used  for 
the  qualitative  analysis  of  rotation  state  errors  but  it  can 
also  be  used  to  carry  out  quantitative  calculations  and  as  a 
guide  for  the  design  of  the  test  system.  The  compressor  rotor 
described  below  is  used  as  an  example: 

Altogether  a  compressor  has  thirty  moving  blades  and 

rotational  speed  n=10700  rotations/min.  Its  flow  field  behind 

the  blade  row  is  as  shown  in  fig.  3a.  T*=238K,  T*=248K,  T  =273K, 

a  d  o 

C2a*93m/sec/  C2b=77m/sec,  1=0.7  and  b=0.3.  The  thermocouple's 

time  constant  T=0.5  sec.  We  know  from  formula  (19)  that  the 
weighted  error  of  the  total  temperature  is: 

&tZ  -  0.7(1  -  0’.7)(284  -  283)(93  7l)/(0.7  X  93  +  0.3  X  77)  —  0.04<C 

A7VAT*-  0.04/(283  -273)-  0.4  %  _ 

(2)  . .  - 

t,  —  40/(10700  X  30)  —4.87  X  lO"1  #(3) 

/*  -  t,  X  5  -  1.87  X-10*4  X  0.3-5.61  X  10-»l>  (4) 

Key:  (1)  Weighted;  (2)  Weighted;  (3)  Seconds; 

(4)  Seconds. 

We  know  from  formula  (23)  that  the  dynamic  response  error  of  the 
total  pressure  is: 

.......  CU _ ... _ _  - . 


....  _  1]  - -0.3=C 

3%.  .  < 

A Tm  -  +  &Tm  -  0.04  -  0.3  —  -0.26<C.  ATa/ATk+—  2.6%. 

(3)  ~  (4)  (5)  (6) 

Key:  (1)  Response;  (2)  Response;  (3)  Rotation; 

(4)  Weighted;  (5)  Response;  (6)  Rotation. 


If  we  only  determine  a  situation  when  b  is  relatively  thin, 
for  example  when  recording  the  design  state  performance,  the 
rotation  state  error  equation  can  be  used  to  find  the  desired 
A  Trotation*^’  The  corresponding  probe  should  have  a  large 


time  constant.  For  the  above  example,  .  ,,  ,=0.04°C,  and 

weighted 

to  cause  ^^rotation*® '  ^rom  formula  (8)  we  know  that 

=  -0.04°C.  By  substituting  this  into  formula  (23), 

rcspons 6 

we  can  derive  T=0.05  sec.  When  a  thermocouple  of  T=0.05  sec. 

is  used  to  determine  the  temperature  rise  of  the  compressor 

rotor  under  this  given  condition,  the  weighted  errors  of  the 

total  pressure  and  the  dynamic  response  errors  are  cancelled  out 

and  A T  .  .  .  =0 . 

rotation 

III.  A  METHOD  FOR  RAISING  MEASUREMENT  ACCURACY  BEHIND  ROTORS 

Rotation  state  errors  and  blockage  errors  are  both  system 
errors.  In  order  to  raise  measurement  accuracy  behind  rotors, 
a  modified  or  controlled  method  can  be  used.  The  main  procedure 
for  reducing  rotation  state  errors  and  blockage  errors  usually  has: 

1.  Corrected  design  probe 

In  order  to  reduce  dynamic  response  errors  and  blockage 
errors,  a  small  probe  should  be  used  as  much  as  possible  under 
conditions  that  fully  ensure  strength  and  rigidity.  Usually, 
the  probe's  cut-off  frequency  should  be  £^=1/2  signal^10  * 

By  using  inertia  compensation,  the  response  of  the  testing 

system  can  be  effectively  raised.  If  the  condition  where  b  is 

relatively  thin  is  only  used  to  determine  the  design  state,  we 

can  rely  on  the  equation  of  rotation  state  error  to  select  a 

suitable  time  constant  to  cause  A T  and  ,  , 

response  weighted 

(or  4Tweighted  and  ^TrotatiQn)  to  cancel  eachother  and  thus 

raise  measurement  accuracy.  The  aerodynamic  blockage  ratio 
should  determine  its  upper  limit  based  on  the  allowable  error 
which  is  generally  less  than  6%.  Using  a  streamline  support, 

"L"  shaped  structure  etc.,  we  can  effectively  reduce  the  aero¬ 
dynamic  ratio.  When  measuring  the  various  parameters  behind  the 
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rotor,  it  is  usually  advantageous  to  use  a  multibranched  small 
single  parameter  probe  of  a  circular  dispersion  point  to  decrease 
the  aerodynamic  blockage  errors.  A  combination  probe  is  suitable 
for  use  where  circular  fields  are  not  homogeneous.  It  is  suit¬ 
able  to  use  a  raked  shaped  combination  probe  when  the  position 
of  the  radial  gradient  of  the  blade  tip,  blade  root  etc.  is 
relatively  large.  If  it  is  only  used  in  places  where  the  average 
mid  radius  and  other  radial  gradients  are  relatively  small,  it  is 
suitable  to  use  a  comb  shaped  combination  probe  to  reduce  aero¬ 
dynamic  blockage. 

2.  Justified  distribution  points 

The  measurement  of  the  section's  mean  parameters  should  be 
realized  by  the  fewest  test  points  arranged  according  to  flow 
field  characteristics.  Because  each  test  point  behind  the  rotor 
is  effected  by  the  mean  working  parameters  of  the  rotor  on  the 
relative  flux,  only  one  or  very  few  test  points  can  be  placed 
along  the  perimeter  behind  the  rotor  when  the  inlet  flow  fields 
are  relatively  homogeneous.  The  axial  position  of  the  test 
points  is  usually  suggested  as  0.2-1. 5b  in  the  literature.  Here, 
we  used  the  relative  measurement  of  chord  length  b  to  show  that 
axial  positioning  is  possibly  because  chord  length  is  a  main 
structural  parameter  of  the  rotor  blade.  However,  chord  length 
and  error  really  have  no  intrinsic  relationship.  From  the 
point  of  view  of  error,  it  is  suitable  to  determine  the  aero¬ 
dynamic  blockage  ratio,  if  conditions  permit,  on  the  basis  of 
permissable  error  and  then  determine  the  axial  position  on  the 
basis  of  the  structural  measurements  of  the  rotor  and  probe. 

3.  By  using  displacement  mechanisms  (rotation,  translation 
and  swing)  and  by  selecting  a  measurement  representative  point, 
aerodynamic  blockage  can  be  effectively  reduced. 
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